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ABSTRACT 



Context. Observations of CO^ suggest column densities on the order 10'^ cm - that can not be reproduced by many chemical models. 
CO^ is more likely to be destroyed than excited in collisions with hydrogen. An anomalous excitation mechanism may thus have 
to be considered when interpreting CO^ observations. Other uncertainties in models are the chemical network, the gas temperature 
or the geometry of the emitting source. Similar is true for other reactive ions that will be observed soon with the Herschel Space 
Obsen'atory. 

Aims. Chemical constraints are explored for observable CO^ abundances. The influence of an anomalous excitation mechanism on 
CO* line intensities is investigated. Model results are compared to observations. 

Methods. Chemical models are used to perform a parameter study of CO* abundances. Line fluxes are calculated for ^(CO*) = 
10'^ cm"^ and different gas densities and temperatures using a non-LTE escape probability method. The chemical formation and 
destruction rates are considered explicitly in the detailed balance equations of the radiative transfer. In addition, the rotational levels 
of CO* are assumed to be excited upon chemical formation according to a. formation temperature. Collisional excitation by atomic 
and molecular hydrogen as well as by electrons is studied for conditions appropriate to dense photon-dominated regions (PDRs) and 
star-forming environments. 

Results. Chemical models are generally able to produce high fractional CO* abundances (x(CO*) ~ 10"'"). In a far-ultraviolet 
(FUV) dominated environment, however, high abundances of CO* are only produced in regions with a Habing field Go ^ 100 and 
r^in <: 600 K, posing a strong constraint on the gas temperature. For gas densities S 10* cm"^^ and temperatures i 600 K, the 
combination of chemical and radiative transfer analysis shows little effect on intensities of CO* lines with upper levels N^p < 3. 
Significantly different line fluxes are calculated with an anomalous excitation mechanism, however, for transitions with higher upper 
levels and densities S 10* cm"^^. The Herschel Space Observatory is able to reveal such effects in the terahertz wavelength regime. 
Ideal objects to observe are protoplanetary disks with densities S 10* cm"^^ . It is finally suggested that the CO* chemistry may be 
well understood and that the abundances observed so far can be explained with a high enough gas temperature and a proper geometry. 

Key words. Astrochemistry - Radiative transfer - ISM: molecules - Submillimeter 



1. Introduction 

Since the first detect ion of the molecul ar ion CO* in the in- 
terstellar medium by iLatter et al] (1 19931) . subsequent observa- 
tions have been carried out successfully toward several pho- 
todissociation regions (PDRs), reflection nebulae and plan- 
etary nebula e (Hoger heiide etalj 119951: IStoerzeret^ 119951: 
Fuente et alj I2003i). pro tostellar envelopes j Ceccarell i et alj 



Stauber et al.ll2007l) and even galaxies (Fuente et al. 20Q3, 
As the = 1 ^ transition is blended by atmospheric 
O2, CO* is most commonly observed in the = 2 — > 1 and 
A = 3 — > 2 rotational transitions. The fine-structure levels 
7 = A + 1 of the lowest few rotational transitions are close in fre- 
quency and they are thus usually observed simultaneously with 
modern spectrometers that offer bandpasses of 1-2 GHz. Higher 
rotational transitions have been obser ved with ISO toward t he 
low-mass protostar IRAS 16293-2422 dCeccareUi et al.lll997h . 

The reactive ion is believed to be sensitive to far-ultraviolet 
(FUV) photons and X-rays and thus to be a tracer for PDRs 
and X-ray dominated regions (XDRs). Chemical models have 
confirmed this by showing that the presence of FUV pho- 
tons or X-ray s enhances the CO* abundance by several orders 
of magnitude (Sternberg & DalgarnollT99l IStauber et al]|2005l: 
[Spaans & Meiierink 2007). The observed CO"^ column densi- 
ties (A(CO*) ~ 10'^ cm"^), however, pose a real challenge to 
PDR and XDR models. They fail to reproduce these abundances 



often by more than an order of magnitude. It should be men- 
tioned though, that 'trad itional' models tread the che mistry ei- 
ther in a plane-parallel (Sternberg & Dalg arndll995) or spher- 
ical ( Stauber et al. 2005) one dimensional geometry. A careful 
treatment of the geometry could therefore be a solution to this 
problem. T he geometry of the emitting source is included in 
the wo rk of iJansen et al.l (1 19951) for the Orion Bar, [Fu ente et jl] 
(I2OO8 ') for the starburst galaxy M 82 and B ruderer et alj ( I2009al) 
for the massive star-forming region AFGL 2591. Their results 
will be discussed in Sect. |2] 



An interesting feature of the observed CO* emission lines is 
that the rotat ional excitation temperatures appear to be as low as 
Tex * 10 K (ILatter et alJll993l:lFuente et alj|2003h . Considering 
also the fact that the timescales for rotational excitation and 
chemical destruction of CO* are similar, suggests that CO* is 
destroyed before its translational motions become thermalized 
(Black 1998). In addition, CO* might be excited upon forma- 
tion. Nascent excitation effects occur when some of the enthalpy 
change in a reaction goes into rotational (or vibrational) excita- 
tion of CO*. On the other hand, if CO* is formed through ioniza- 
tion of CO by energetic photons or particles, it will tend to inherit 
the rotational excitation of its parent because the relatively heavy 
nuclei can not respond during the rapid electronic ionization pro- 
cess. It is therefore not clear how the rotational CO* levels are 
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excited and what influence an anomalous excitation mechanism 
has on the hne fluxes. 

Motivated by the shortcomings of traditional PDR and XDR 
models and the considerations in the previous paragraph, the aim 
of this paper is to study the excitation of CO^ by combining 
chemical models and radiative transfer analysis. Uncertainties in 
chemical models are the gas temperature and the chemical net- 
work. Thus, we first explore the chemistry of CO^ in a general 
parameter study to constrain the physical gas conditions required 
to produce large fractional abundances (Sect.|2]i. We then make 
use of RADEX, a compute r program for fast non- LTE analysis 
of interstellar line spectra ( van der Tak et a l. 2007), to calculate 
the rotational excitation of CO"^, including the formation and de- 
struction rates explicitly. The method is described in Sect. 13.11 
and l3.2l Calculated line fluxes for selected frequencies are pre- 
sented and discussed in Sect. 13.31 The conclusions of this study 
are drawn in Sect. |4] 



2. Chemistry 

The chemical network for CO ^ has been d iscussed in the 
past by many authors (e.g., | Ster nberg_& Dalgarn o, |l995t 
ISavage & Ziurvsl 12004: Stauber et al. 2004, 2005, |2^ and 
most recently by iBruderer et al., t2009a) . Therefore, we give 
only a short summary of the most important reactions. 

Under the influence of FUV photons, CO^ is efficiently 
formed in reactions of with OH. stems from photodis- 
sociation and subsequent ionization of CO and C, respectively. 
OH has high abundances at temperatures above ^ 300 K, where 
large amounts of oxygen are driven into OH due to the disso- 
ciation of gaseous water. For low FUV fluxes, other formation 
routes are the reactions of C^ with CO2 and O2. CO^ rapidly re- 
acts with H2 and H to form HCO+, HOC+ and CO, respectively. 
Another fast destruction mechanism is the dissociative recombi- 
nation with electrons. 



2.1. Parameter study 

The strong dependence on the OH abundance and the fact, that 
CO^ is quickly destroyed by hydrogen makes its abundance very 
sensitive to the gas temperatur e and density. To show this, we 
have used the chemical grid of Bruderer et al ] (I2009 b'> to calcu- 
late the fractional abundances of CO^ for different gas tempera- 
tures Tkin, gas densities n(H2) and FUV fluxes Gp (in units of the 
Habing field). The models of iPellegrini et alJ (l2009i) for Orion 
Bar showed that CO^ is also sensitive to cosmic rays. We do not 
further investigate their influence but include a generic cosmic 
ray ionization rate of 5.6 x 10"'^ s"' . It should be noted that the 
chemical equilibrium of CO^ for the densities and temperatures 
mentioned above is reached very quickly. Our time-dependent 
chemical models indicate that the CO^ abundance remains con- 
stant after x 1000 years. 

Figure [T] shows the results for gas densities n(H2) - 10"*- 
10^ cm^^ and an optical depth Ay = 0.5. At this optical depth, 
CO'^ is b elieved to reach its maximum abundance in dense PDRs 
( I Sternbe rg & Dalgarno 1995). The observed column densities 
are between ~ 10"-10'^ cm"^ (Table©. Assuming all C0+ to 
be produced at Ay = 0.5, these column densities coiTespond to 
fractional abundances of x(CO^) ~ lO '^-lO"^. According to 
Fig.n x(CO+) * lO ^* is produced for n(H2) = 10^ cm^^ wifli 
Tidn ^ 1000 K and Go ~ 10^^-10^ for n(H2) = 10^ cm^^ wifli 
Tidn ^ 1600 K and Go ~ 10^-10^ and for n(H2) = 10*^ cm^^ 
with Tkin ^ 1000 K and Go ^ 5 x 10"*. A fractional abundance of 



10"'" can be reached for n(H2) = 10"* cm^^ with r^n ^ 600 K 
and Go ~ 10^-10^; for «(H2) = 10^ cm^^ with Tkin ^ 600 K and 
Go ^ 10^-5 X 10^ and for n(H2) = 10^ cm^^ with T^in ^ 600 K 
and Go S 10^. For a typical FUV flux Go ~ 10"* in PDRs, abun- 
dant CO^ (x(CO^) = 10 can therefore only be expected in 
regions with gas temperatures Tidn ^ 600 K. 



2.2. Comparison with observation 

Although the chemical model results presented in Sect. 12.11 are 
completely independent of geometry, a simple comparison with 
observations can be made by assuming the clouds to be large 
enough in size to provide the observed column densities for 
Av = 0.5. The assumption is reasonable since the region needs 
to be only 1000-10000 AU in size for «(H2) = 10^ cm^^ and 
A^(CO^) = 10"-10'^cm"^ whereas most PDRs are larger (de- 
pending on the viewing angle of course). Table [1] lists CO^ col- 
umn densities observed toward well known PDRs. Also shown 
is the ratio of the column density to a hydrogen column density 
at Ay = 0.5 (A^(H2) ~ 2 x 10^' cm"^). This ratio allows direct 
comparison with the chemical models. It is apparent that the ob- 
servations are in fair agreement with the model results presented 
in Fig.[T] For example, M17SW and Orion Bar have a fractional 
abundance of a 10"'' with N(CO+) ~ 10'^ cm^^ at Ay = 0.5 
in regions with n(H2) ~ 10^-10*' cm"^ and Go = 5 x lO"*. Such 
an abundance is consistent with our models for Tki,, <: 1000 K. 
NGC 7023 has a somewhat smaller gas density and Go = 2x10"'. 
The observed abundance is comparable to our models for Tkin ^ 
600 K. Mon R2 and G29.96-0.02 on the other hand are believed 
to have gas densities n(H2) ~ 10^ cm"^ and Go ~ 10^. The col- 
umn densities of these dense objects are in agreement with our 
models for ^ 800 K. 



2.3. Discussion 

The parameter study presented above is a simplification of the 
physics and chemistry in order to study the CO^ abundance for 
a few eff'ects such as FUV flux, temperature and density. In re- 
ality, the conditions for the chemistr y is expected to vary from 
one source to another. For example. Pellegrini et al.l ( 12009 ) mod- 
eled the Orion Bar P DR using the spectral synthesis code Cloudy 
dFerland et al.lll998 ) to derive the physical conditions across the 
Bar. They assume the cosmic rays to be trapped in the cloud by 
a tangled magnetic field which eventually yields a highly en- 
hanced cosmic ray ionization rate. This allows them to success- 
fully reproduce the observed CO^ column densities. However, 
the simple comparison in the previous section between models 
and observations shows good agreement if a certain temperature 
is assumed. The problem of traditional PDR models to produce 
high CO^ abundances may thus not lie in the chemistry itself but 
in the uncertainty of the gas temperature at Ay ^ 0.5-1 which 
needs to be ^ 600 K. The gas temperature is usually calculated 
self-consistently in these models by solving the thermal balance. 
The various PDR codes, however, show diff'erences on the order 
of a magnitude for Ay = 0.5-1, ranging fro m a few hundred 
to a few thousand Kelvin dRollig et al.ll2007b . This uncertainty 
results eventually in different CO^ abundances. CO^ may thus 
serve as a tracer for high gas temperatures at low Ay. 

The main uncertainty in the chemical models, on the other 
hand, are the reaction rates. For example, the rate coefficient for 
the C^ + OH reaction in our models is 7.7 x 10"'*' cm"-' s"', 
independent of temperature (UMIST, Woodall et al. 2007 ). The 
temperature dependent rates published bv .Dubernet et al.l (1 19921) 
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Fig. 1. Contour levels for the CO^ fractional abundances for hydrogen densities between IC* cm ^ and 10^ cm ^. The assumed 
optical depth is Ay - 0.5. The contour levels are equally spaced from 10"''* to 10"'" and extended by 5 x 10"'" and lO"**. 



Table 1. CO^ column densities derived from observations to- 
ward PDRs 



Object 


"(Hz) 


Go 


A' 


N/2 X 10^' 




[cm"3] 


[Habing] 


[10'- cm"2] 


[10-'] 


M17SW;i 


5 X lO" 


5x 10* 


1-1.8 


0.5-0.9 


Orion Bar|^ 


2.5 X 10^ 


5 X 10* 


1-2.7 


0.5-1.35 


8140;. 


10= 


150 


0.03 


0.015 


NGC70231 


lO^-lO^ 


2x 10^ 


0.3 


0.15 


Mon R2^ ~ 


1.5 X 10'' 


4.9 X 10^ 


0.53 


0.27 


G29.96-4).02'' 


7 X 10^ 


1.5 X 10^ 


0.47 


0.24 



" iLatter et al.' ('l993'):'Stoerzer et al.' l!l995h 

* ^toerzer et al. (1995); Fuente et al. ( 2 003b;ISavage & Ziurvsl ([2004h 
" IPark & Minhl (119951) : iStoerzer et al] ( Il995h : ISavage & Ziurvsl 

( 120041) 

IFuente & Martin-Pintadol jl997b ; lFuente et al.l j2003l) 
" FrIzzo et al] ( f2003l) 



and ItYo3 (Il996h on the other hand, are above 10 cm s . 
However, the branching ratio is not clear as the + OH — > 
CO + H^ reaction is somewhat more exoergic than the CO^ + 
H channel. It is interesting, that the former reaction, producing 
H"'", may partly recycle OH, because H^ has a high probability 
of charge transfer with O in neutral gas at the temperatures of 
interest and the resulting O^ will lead b ack to OH and H7O as 



long as the H2/H ratio is not too low (e.g . ,ISternberg & Dalgarnol 
[1995; Stauber et al. 2006). Fuente et al] (l2008h were able to re- 
produce the observed CO^ abundances with a reaction rate of 
2.9 X 10"'' cm"^* s"' . Assuming that the abundance roughly scales 
with this reaction rate, the modeled CO^ abundances may well 
be 3^ times higher. 

A strong support for the argument that the CO^ chemistry 
may actu ally be well underst ood comes also from recent model 
results by 'Fuente etaTI (l2008l) and lBruderer et all ( l2009al) where 
the geometry of the observed object has been taken into account. 
iFuente et al.l (2008) modeled M 82 by assuming a plane-parallel 
cloud illuminated by FUV photons from two sides. With this 
more realistic view of the interstellar medium in M 82, they 
obtained ^ 5x larger CO^ column densities in good agree- 
ment with the observations and increased the go odness of fit 
for other observed species. Br uderer et al. i (f20()9ah constructed 
a two dimensional model of the massive star-forming region 



AFGL 2591 to include FUV irradiated outflow walls. They 
found that the FUV fluxes and gas temperatures along the walls 
are high enough to produce large amounts of CO^ . Their models 
were able to reproduce the observed CO^ abundance within a 
factor of ^ 2, whereas the one dimensional spherical models of 
IStauber et all (|2004 l2005h underestimated the CO+ abundance 
by several orders of magnitude. A reason for this is that a suitable 
geometry increases the Ay ~ 0.5 area along the line of sight. 

X-r ays are able to enha nce the CO^ abundance to similar 
values. IStauber et all ( l2007h and lBruderer et al] ( l2009ah showed 
that X-ray fluxes Fx ^0.1 erg s"' lead to fractional abundances 
x(CO^) ~ 10""'-10"^. However, since most observations of 
CO^ are believed to trace FUV photon s rather than X-rays (e.g., 
IStauber et al]|2007t iFuente et al]|2008i) . the chemistry of CO+ 
in XDRs is not explored and we refer to the papers mentioned 
above for further information. EUV photons with energies above 
14 eV might be an important source for CO^ at the boundaries 
of photoionized nebula e. This is not included in the models of 
iBruderer et al] (l2009bh though. EUV photons are quickly ab- 
sorbed by the gas and the effect on the total CO^ abundance 
is expected to be small. However, EUV chemistry would only 
increase the CO^ abundance and thus improve the goodness of 
fit of chemical models when compared to observations. 

Although a careful treatment of the gas temperature and the 
geometry of the object is crucial to interpret CO^ observations 
properly, the question remains what the influence of the excita- 
tion mechanism is on emission lines. This will be studied in the 
next section. 



3. Radiative transfer analysis 

3.1. Method 

Since CO^ is destroyed by hydrogen and electrons on a rela- 
tively short timescale, the chemical formation and destruction 
rates need to be considered when calc ulating the statistical equi- 
librium (e.g. . Ivan der Tak et a D2007h : 



drii 
dt 



(1) 



where is the number of levels considered in the model, n, 
[cm"-'] is the level population of level /, Fj and D, [cm""* s"'] 
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are the chemical formation and destruction rates for level /, re- 
spectively, and Pjj [s^^] is given by 



Aij + BijJ, + Cij (Ei > Ej) 



(Ei < Ej) 



(2) 



Aij and Bij are the Einstein coefficients for spontaneous and in- 
duced emission, Cij are the rate coefficients for collisions, Jy is 
the specifice intensity integrated over the line profile and aver- 
aged over all directions and £, is the energy of level /. 

Following the results of Sect.|2l we assume that the chemical 
rates are in equilibrium {F = 2, - 2, D, = D). It is further 
assumed that CO^ is formed in an excited state, where the levels 
are populated according to a formation temperature Tform and 
that all levels are destroyed with equal probability. Hence, the 
chemical formation and destruction rate for level / are 



-E,/kT!, 



Ft = 



ljj=l 1 

D 

Di = - , 
N 



^//^^form 



(3) 



(4) 



where gi is the statistical weight of level /. Equation [T] and the 
equilibrium between formation and destruction (F = D) yield 



= Xi«>^./7-«,£A7 + Z) 



1 



(5) 



It should be emphasized that the formation temperature is an ar- 
tifice that allows to describe the effect of the nascent population 
distribution by a single parameter. The formation temperature 
should not be conflated with the physical temperature nor with 
the excitation temperatureB 

For a consistent solution of these equations, the condition 
«/ = n(CO^) needs to be introduced. It is convenient to 
scale Eq. |5] with l/n(CO^) since the equations become inde- 
pendent of n(CO^). This scaling does not change the excitation 
temperature and hence the line intensity obtained from the level 
populations. It is further justified by the assumption of F - D 
which ensures that the chemical rates are proportional to n(CO^) 
(see also Eq. |7]i. However, the system of equations can now be 
solved for the fractional abundances of levels (jc, = «,/n(CO^)): 



T.%gje-^'l"<^ 



N 



, (6) 



with the normaUzed destruction rate 



(7) 

where n(i) is the abundance [cm of a molecule i and kiiTi^^) 
is a rate coefficient [cm^ s"^] depending on the kinetic temper- 
ature. Due to the high abundance of electrons, atomic hydrogen 
and H2, we can savely take the sum only over these species to 
calculate the destruction rate: 

D ~ 1.5x10A(H2)h-7.5x10"'VH)+2.0x10"^(:^^)«(O ■ 

\ 300K/ 

(8) 



' The excitation temperature of a transition in the radiative transfer 
model is defined by the ratio of the upper level population to the lower 



The rate coeffi cients are taken from the UMIST database for 
astrochemistry dWoodall et al.ll2007h . According to UMIST, the 
accuracy of these rates is between 25% and 50%. 

To calculate the intensities for a range of rotational transi- 
tions, we mak e use of the publicly av ailable radiative transfer 
code RADEX jvan der Tak et alj2007l) . RADEX allows to com- 
pute quickly the intensities for a grid of defined gas properties in 
a uniform medium, based on statistical equilibrium calculations. 
The program includes radiation from background sources and 
treats optical depth effects with an escape probability method. 



one, that is by x,+i/x; = g;+i/^,e 



-hv/Ta 



For all calculations, we assume a CO^ column density of 
A^(CO+) = 10'^ cm-^, a line width of AV = 1 km s^' and a 
2.73 K blackbody background. With these values, the lines dis- 
cussed in the following sections are optically thin. 



3.2. Molecular data 



The energy levels and transition frequenc ies of CO^ are taken 
from the JPL database (iPickett et al .1119981) . The transition prob- 
abilities are calculated for a dipole moment fi - 2.63 D 
(ICheng et al. 2007). Since no published excitation rates are 
available for_CO^-H2 collisions, we apply the HCO^-H2 rates 
of iFlowed (IT999.) . These rates were calculated for temperatures 
between 10 K and 4 00 K and rotational levels up to 7 = 20. 
ISchoier et al.l (12005 *) extrapolated this set of coefficients to in- 
clude energy levels up to 7 = 30 and temperatures up to 2000 K. 
We adapt the same rates for the fine-structure levels. In contrast 
to C0+ (X HCO+ {X '2;+) has a 49% higher dipole mo- 
ment and a 24% smaller rotational constant. In addition, the H2 
excitation rates are likely to be different due to the nonzero spin 
of CO^. How ever, the error may not be bigger than factors 2 - 3 
(see e.g.. Bla ck & van Dishoeck|[T99l1; ISchoier et al.ll2005h . For 
rotational transitions induced by ato mic hydrogen, we use the H- 
CO^ deexcitation rates calculated bv lAndersson et all (l2008l) for 
y < 8. The rates have been extrapola ted to levels up to / = 30 
following the method as described bv lSchoier et al.l (l2005l) . The 
rate coeffi cients for electron deexcita tion are calculated using 
Eq. 2.9 in [Dickinson & Flqw ei* (1981 1). Since the fine- structure 
is not treated properly bv [Dickinson & Floweij (Il98ll) . we as- 
sume the sa me rates for the corresp onding levels. It should be 
noted, that iFaure & Tennvsonl(l200ll) have investigated the elec- 
tron impact excitation of CO^ in more detail. Although they 
neglect the fine-structure too, it is interesting that the rates for 
AA^ > 1 transitions are significant, something that is not antic- 
ipated by Dickinson & Flower (1981). However, since the for- 
mula of Dickinson & Flower. (. 198 1) allows to calculate the rates 
also for high rotational levels, and since electron impact excita- 
tion is of minor importance for our purpose, we adopt the treat- 
ment of Dickinson & Flower ( 198l])- 

The critical densities (ncrit - Aui/YjKui) of our applied 
molecular data for Tkin - 500 K are presented in Table |2] 
The lines in the table are chosen due to their observability with 
ground based telescopes. It can be seen from the critical densi- 
ties that collisions with atomic hydrogen become important for 
H/H2 ratios of ^ 1 . Electron excitation will be important for an 
ionization fraction x(e^) k. 10"^. The critical electron densities 
are in reaso nable agreement (differenc e is $ 20%) with those 
published by lFaure & TennysonI (l200ll) for Tkin = 500 K. 
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Table 2. Critical densities for r^n = 500 K. 



Species 


A' = 2 ^ 1 
«c,it [cm-3] 


N = 3^2 
n,,ii [cm-3] 


A^ = 4^ 3 
"crit [cm-^^] 


A' = 7 ^ 6 
'icit [cm-3] 


Hj 


4.2 X 10' 


1.1 X 10" 


2.1 X 10" 


8.4 X 10" 


H 


6.5 X 10' 


1.6 X 10*^ 


3.0 X 10" 


1.4 X 10' 


e" 


3.9 X 10- 


1.4 X 10^ 


3.4 X 10' 


1.9 X 10^ 



3.3. Results and discussion 
3.3.1 . Line fluxes 

The integrated line intensities are calculated for a grid of gas 
densities (n(H2) = 10'*-10^ cm"-'), gas temperatures (Tidn - 10- 
2000 K) and formation temperatures (rform = 10-3000 K). 
Although the line fluxes will also depend on the chemical forma- 
tion and destruction rates, we treat them as as being well known 
(Eq. and do not vary the chemical rate coefficients, for in- 
stance. Also calculated are the fluxes for £) = 0, that is the radia- 
tive transfer models without taking the chemical processes into 
account. To study the influence of atomic hydrogen and electron 
excitation, a model was set up with H/H2 = 1 and x(e~) = 10"^. 
These are typical values for PDRs with n(H2) 10^ cm"^ and 
Go ~ 10^. It is found that the results for the = 2 ^ 1 lines 
differ 50% at most from those where only H2 was considered. 
For simplicity, atomic hydrogen and electron excitation is thus 
neglected in the following paragraphs. We have also carried out 
calculations considering the dust temperature. Models were run 
with the dust temperature varying from 20 K to 1500 K with 
a gas density of 10^ cm"-' and a total H2 column density of 
10^^ cm"^. The gas temperature was assumed to be the same 
as the dust temperature. It was found that the influence of dust 
on the CO^ line fluxes studied here can be neglected under these 
conditions. Since most observed lines are those with = 2 ^ 1 
and A^ = 3 — > 2, we will concentrate on these transitions. 

Figure |2] shows the integrated line intensities (K km s"') for 
the 2|^ ^ l| transition at 236.06 GHz for different formation 
temperatures. Also shown is the flux assuming LTE conditions. 
The LTE flux was calciilated using the analytical expression 
given in 1 Stauber et al.l (12007 ). As expected, the LTE fluxes are 
only matched by models with 2) = (no formation temperature 
assumed) for n(H2) ^ 10^ cm"^, in other words for densities 
where collisions dominate and the levels become thermalized. 
The assumption of a formation temperature Tform ~ 500 K on 
the other hand implicates that the levels do not become thermal- 
ized (see plot for n(H2) = 10^ cm"^ in Fig.|2]i. 

The results for the different formation temperatures vary 
only ^ 20-30% for low densities (n(H2) $5x10^ cm"^) and 
temperatures where CO^ has high abundances (r^in ^ 600 K; 
shaded region in the figure). At higher densities, the models in- 
cluding the chemical formation and destruction rates show fac- 
tors 2-7 higher line fluxes. Models with formation temperatures 
between 500 K and 1000 K have similar fluxes like those with 
2) = 0. This is also true for models with formation temperatures 
exceeding 1000 K for n(H2) 5 10' cm"^. At higher densities, 
the upper levels of the A^ = 2 — > 1 lines are not or only weakly 
populated for formation temperatures exceeding 1000 K. 

Figure [3] shows the results for the 3^ ~* 2| transition at 
354.01 GHz. The fluxes for the different formation tempera- 
tures vary only w 2Q-A0% for densities ^10^ cm""*. At higher 
densities, the fluxes can be up to 3 times higher, depending on 
the formation temperature. Models with formation temperatures 
^ 500 K have similar line fluxes as models where 2) = 0. It 
can be concluded that the influence of an anomalous excitation 



temperature on the line flux is small for the CO^ transitions and 
regions observed so far (n(H2) ^ 10^ cm"-'). 

Bigger differences may be expected for transitions with up- 
per energy levels comparable to the gas temperature where CO^ 
is most abundant. The high spectral resolution instrument HIFI 
on board the Herschel Space Observatory covers the frequency 
range 480-1910 GHz and thus the upper CO+ levels A^up = 5- 
16 with upper energy levels between 85 K and 769 K. Figure |4] 
shows the line fluxes for transitions N —> N-l of the first 16 up- 
per levels for n(H2) =10** cm""* and Tkin - 600 K. For the sake 
of clarity, only the stronger one of the two fine-structure lines are 
presented. The density of 10** cm"^ represents the critical den- 
sity of the higher A^up-levels. It can be seen that the various for- 
mation temperatures result in significantly different line fluxes. 
For example, transitions with A^up > 14 have an order of magni- 
tude lower line fluxes when the formation temperature is much 
lower than the coiTesponding upper energy level (rfonn ~ 10 K). 
Further calculations for n(H2) = 10^ cm"^ show that this is al- 
ready the case for A^up > 9. Calculations for n(H2) = 10"* cm"^, 
on the other hand, indicate that transitions with A'up ^ 7 have 
line fluxes below 10 mK and may there fore not be obs e rvable 
for A^(CO+) = 10'- cm ^ Interestinglv. iCeccareUi et al.l (1 19971) 
detected far-infrared CO^ lines up to A^ - 21 towards the low- 
mass star-forming region IRAS 16293-2422. Clearly, the for- 
mation temperature of these CO^ levels is either high or there is 
no nascent excitation due to direct ionization of CO. This might 
indicate the presence of X-rays, enhanced cosmic rays or EUV 
photons. Nevertheless, dense PDRs and protoplanetary disk at- 
mospheres with high gas densities (^ 10^ cm"') and high tem- 
peratures (<; 600 K) are ideal testbeds to study the excitation of 
CO^ in the terahertz wavelength regime with Herschel. 

Similar might be true for other molecular ions that will 
be observed with Herschel. CH+, for example, is als o mor e 
likely to be destroyed than excited by hydrogen (Black I T998I) . 
The chemistry for CH+, however, is more complex. The main 
source of CH^ is the endoergic reaction C^ + H2 ^ CH^ + H. 
Besides high gas temperatures, vibrationally excited hydrogen 
could carry a significa nt part of the ener gy needed to activate 
the reaction (Sternberg & Dalgarnolll995h . This introduces an- 
other uncertainty to chemical models. In addition, CHLis found 
to be sensitive to the thermal emission by dust (lBlacklll998l) . 
Nevertheless, care should be taken when analyzing observations 
of such reactive ions. They do not necessarily reflect the column 
density and hence the molecular abundance when they are ex- 
cited upon formation. 

3.3.2. Excitation temperatures 

It is interesting to see that the calculated line fluxes for gas den- 
sities and gas temperatures of typical dense PDRs («(H) of a few 
xlO*" cm""* and Tkin ^ 600 K) are comparable to the maximal 
LTE line fluxes (see LTE fluxes at T]^i„ = 18 K and Tex = 35 K 
in Fig. 12] and Fig. [3] respectively). For these conditions, the as- 
sumption of a low excitation temperature therefore yields ap- 
proximately the same flux as detailed non-LTE calculations. 

It has already been mentioned, that the observed CO^ 
lines sh ow rather low (Tex ~ 10 K ) excitation temperatures 
(iLatter et all 1 9931; iFuente et al.l2003h . This agrees well with our 
models. The calculated non-LTE excitation temperatures for a 
CO^ column density of 10'^ cm"^ and n(H2) = 10^ cm""* are 
between w 15 K and » 70 K (Fig. |5j for gas temperatures be- 
tween 600 K and 2000 K. The excitation temperatures for the 
N - 3 —> 2 line are equally low (Tex ~ 10-25 K). They are 
clearly below the gas temperature and thus far from LTE. 
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Fig. 2. Integrated line intensities for the 2| — » 1 1 transition at 236 GHz. The soHd Hne corresponds to the radiative transfer models 
with D = (no formation temperature is assumed). The vertical line at Tkin = 18 K indicates the position of the peak LTE flux. The 
shaded region marks the gas temperatures where CO^ is chemically most abundant. 



The assumption of a formation temperature with rfo„n — 
10 K leads to an excitation temperature that matches the ob- 
served one best. T his s uggests that the CO ^ lines observed by 
iLatter et all (Il993h and iFuente et all ( |2003|) were excited upon 
formation with an excitation energy corresponding to a low for- 
mation temperature. 

3.3.3. Line ratios 

Line ratios are useful to estimate physical gas conditions. 
Figure |6] shows the contour lines for different densities for the 
J - 2 J ll/^l — » 2| flux ratio as a function of the gas 
temperature and formation temperature. The contour lines for 
the flux ratio depending only on gas temperature and density 
(D - 0) are given in Fig.|2l It is seen in general that the depen- 
dence on the formation temperature is strongest for low gas tem- 
peratures. The results for the models without taking the chemical 
formation into account are similar to those with high formation 
temperatures. 

Comparison of Figs. |6] and [T] reveals that ratios 7 = 2| — > 

l|/3| ^ 2f > 3 are only possible for «(H2) ~ 10"* cm-^ 
Tkin < 1000 K and formation temperatures between 10 K and 
100 K. Such high ratios are thus indicative for an anomalous 
excitation mechanism with low formation temperatures. 

The flux ratios for the two = 2 ^ 1 and = 3 — > 2 
fine structure lines are x 0.5 and x 0.6-0.7, respectively. Their 
dependence on gas temperature and density is weak as expected 
in the optically thin limit. In addition, the ratios differ at most 
20% for various formation temperatures and are thus not suitable 
to estimate the gas conditions or formation temperature. 



3.3.4. Comparison witli observations 

In general, the derived column densities from the observed line 
fluxes are consistent with ou r model resul t s pres ented in the pre- 
vious section. For example. IFuente et al.l (|2003|) obtained a col- 
umn density of A^(CO^) - 10^^ cm~^ for an observed flux of 
0.63 K km s ' toward Orion Bar. This flux is comparable to the 
radiative transfer models for n{H2) ~ 10^ cm"^, Tidn ~ 600 K 
and Tform ~ 10 K. 

However, when comparing observations with the homoge- 
neous radiative transfer models described in the previous sec- 
tion, four parameters need to be fitted: the column density, 
the formation temperature and the gas density and temperature. 
Therefore, detections of at least four different lines are necessary. 
Unfortunately, more than two observed CO^ transitions are re- 
ported only for Ori on Bar (Hogerheijde et al. 1995) and M17SW 
^Latter et al.lll993l) . Nevertheless, the CO^ fluxes observed to- 
w ard these two regions are modeled using a test as described 
in iHogerheiide et akl d 19951) . The set of parameters which min- 
imizes is then considered as best fit to the data. Since the 
critical densities are similar for the lines that are accessible by 
ground based telescopes, CO^ is not the ideal molecule to trace 
the gas density. The density is thus adapted as described in the 
following paragraphs. 

The molecular cloud Ml 7SW lies to the southwest within the 
H II region of the Omega Nebula. T he PDR is separating t he ion- 
ized gas from the molecular cloud. iMeixner et al.l (1 19921) found 
that the PDR consists of clumps with n(H2) = 2.5 x 10^ cm"^ and 
T X 1000 K and a lower density core surrounding the clumps 
with n(H2) = 1500 cm"^ and T x 200 K. The same authors es- 
timated a FUV field of Go ~ 5 x lO''. In a homogeneous model. 
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Fig. 3. Integrated line intensities for the 3 5 — » 2| transition at 354 GHz. The soHd Hne corresponds to the radiative transfer models 
with D = (no formation temperature is assumed). The vertical line at Tkin = 35 K indicates the position of the peak LTE flux. The 
shaded region marks the gas temperatures where CO* is chemically most abundant. 



the density was found to be 1.5 x 10"* cm""* and the tempera- 
ture 300 K. Thus, for the radiative transfer models, a density of 
«(H2) = 1.5 X 10"* cm""* is assumed with a H/H2 ratio of 1 and 
an electron fraction x(e ) - 10 The error in the observations 
for the test is taken to be 30% of the observed flux for the 
/ = 2| ^ 1 1 lines and 50% for the others. The fits with;i'^ < 1 
have column densities A^(CO*) = 1.5 x 10'^^.5 x lO'^ cm"^ 
and gas temperatures <: 100 K in good agreement with the val- 
ues presented in Table [T] Unfortunately, no conclusions can be 
drawn regarding the formation temperature. However, the flux 
ratio of the observed J = 2 | — » l| and J - 31 ^ 2| lines 
including the error is » 5-8 dLatter et al. Further calcu- 

lations show that this indicates gas temperatures $ 800 K and a 
formation temperature between 10-70 K. 

The Orion Bar is a dense molecular ridge within the Orion 
Molecular Cloud illuminated by the Trape zium stars. From 
millim eter and submillimeter observations, 'Hoge rheiide et alj 
(Il995h find that the Orion Bar is best described by a clumpy 
medium where x 10% of the material may be in clumps with 
n(H2) ~ 10^ cm"^ and » 90% in a homogeneous interclump 
medium with n(H2) ~ 3 x 10"^ cm~^. The FU V field is esti- 
mated to be Go = 3 X 10"^ dJansenet all 1 1995 '). To model the 
observed CO* emission, we use a weighted average density of 
1.3 X 10^ cm-^ with a H/H2 ratio of 1 and xCe") = 10 The 
line widths and fluxes for the two = 2 — > 1 and N - 3 
2 transitions, respective ly, are taken from the observations of 
iHogerheiide etall (11995'). The error for the test is taken to be 
30% for the = 2 ^ 1 transitions and 50% for the = 3 ^ 2 
lines since these pr ofiles are dominate d by in strumental broad- 
ening according to IHogerheiide et al.l (Il995h . Best fit models 
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< 0.8) show column densities A^(CO*) = 5.5 x 10' 



10" cm ^ and gas temperatures r^in ^ 700 K. The mod- 



els, however, are not conclusive regarding the formation tem- 
perature. Nevertheless, the fits are in good agreement with the 
chemical model results in Sect. |2] and with previously reported 
column densities (Table [T]i. 

The low observed excitation temperatures and the high J = 
2|_-» l|/3f ^ 2| flux ratio ( Latter et al.lll993l; iFuente et al.l 
I2003h indicate that CO* may be formed rotationally excited with 
a low formation temperature (rform ~ 70 K). Due to the as- 
sumptions made in the models (HCO*-H2 rates, homogeneous 
medium) we consider this as indication rather than proof though. 
One should also bear in mind that the chemistry and excitation 
conditions may vary from source to source. If CO* is mainly 
produced by ionizations of CO, no nascent excitation effects will 
occu r as stated in Sect. [T ] This is the cas e by the presence of X - 
rays dStauber et al.ll2005h or cosmic rays dPellegrini et alj|2009l) . 
UV photons will have this effect only at the boundary of a PDR 
and is therefore expected to be small. Certainly, more observa- 
tions of different transitions are needed toward PDRs to further 
investigate this problem. In addition, multi-dimensional radia- 
tive transfer models in combination with chemical models are 
necessary to deal with the complex geometry of the emitting ob- 
jects (Bruderer et al., in preparation). 



4. Conclusion 

A grid of time-dependent chemical models has been calculated 
to constrain the physical parameters where CO* is most abun- 
dant. Line fluxes have been computed for A^(CO*) 
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by tfie use of a non-LTE radiative transfer code to study tfie ef- 
fects of an anomalous excitation mecfianism. Tfie following list 
summarizes tfie main conclusions to be drawn from tfiis study: 

1. Higfi fractional abundances (x(CO*) ^ 10 '° for Ay - 
0.5) are only reacfied for gas temperatures k. 600 K. A sim- 
ple comparison between cfiemical models and observations 
sfiows good agreement if a certain gas temperature is as- 
sumed. Tfiis suggests tfiat tfie CO* cfiemistry is well under- 
stood and tfiat tfie molecular ion serves as a tracer for fiot gas 
in regions witfi Ay ~ 0.5 (Sect. lTTT i. 

2. Tfie model results in Sect. 13.3.11 indicate tfiat tfie influ- 
ence of an anomalous excitation mecfianism is small on 

= 2 — > 1 and N = 3 —> 2 line fluxes in regions wfiere CO* 
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fias been observed tfius far (PDRs witfi n(H2) ~ 10^ cm"^). 
Tfie calculated line fluxes are in good agreement witfi obser- 
vations and confirm tfie cfiemical model results (Sect. [334l i. 
3. Formation temperatures in tfie range of 10-1000 K sfiow 
significantly different results for transitions witfi A^up ^ 9 
for «(H2)_= lO*" cm-3 and A^up > 14 for «(H2) = lO** cm^^ 
(Sect. |3.3Tn i. Tfie Herscfiel Space Observatory covers tfie up- 
per CO* levels A^up - 5-16 and is tfius ideally suited to study 
tfie rotational excitation of CO* and related molecular ions. 
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4. Rotational CO* levels with temperatures much higher than 
the formation temperature {Eap/k » rform) will be scarcely 
populated when excited upon formation. Transitions from 
these levels will therefore have low line fluxes and may not 
be detectable (Sect. 13.3.11 ). 

5. The low excitation temperatures which are observed 
for CO* are consistent with the radiative transfer models 
(Sect. UXSl i. Comparison with observations suggest forma- 
tion temperatures of only 10 K. 

6. The N - 3 —> 2/2 —> 1 line ratios are found to be enhanced 
in Sect. l3.3.3l for low formation temperatures (rform 5 70 K). 
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